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We study the heavy-heavy-light quark (QQq) potential in SU(3) quenched lattice QCD, and 
discuss one of the roles of the finite-mass valence quark in the inter-quark potential. Monte Carlo 
simulations are performed with the standard gauge action on the 16 4 lattice at /3 = 6.0 and the 
0(a)-improved Wilson fermion action at four hopping parameters. For statistical improvement, the 
gauge configuration is fixed with the Coulomb gauge. We calculate the potential energy of QQq 
systems as a function of the inter-heavy-quark distance R in the range of R < 0.8 fm. The QQq 
potential is well described with a Coulomb plus linear potential, and the effective string tension 
between the two heavy quarks is significantly smaller than the string tension a ~ 0.89 GeV/fm. 
It would generally hold that the effect of the finite-mass valence quark reduces the inter-two-quark 
confinement force in baryons. 

PACS numbers: 11.15.Ha, 12.38.Gc, 14.20.Lq, 14.20.Mr 



I. INTRODUCTION 

In hadron physics, the inter-quark interaction is one of 
the fundamental and essential properties. In particular, 
the quark confinement is not only an important property 
of hadrons but also an important problem of the modern 
physics. The non-Abelian and strong coupling nature of 
QCD makes it difficult to treat the inter-quark interac- 
tion and many other nonperturbative phenomena ana- 
lytically. Lattice QCD is the first-principle calculation 
based on the QCD Lagrangian, and is one of the most 
useful approaches for such nonperturbative phenomena 

The quark confinement in hadrons is well described 
by the picture of the gluonic "flux tube" or "string" [3|. 
This means that the confinement potential is a linear 
function of the flux-tube length. For example, the quark- 
antiquark (QQ) potential is written as the one-gluon- 
exchange Coulomb potential plus the linear confinement 
potential, 



Vqq(R) = <tqqR- 



R 



Cnn- 



(1) 



which is called the Cornell potential Here R is the 
distance between the quark and the antiquark, and it is 
equal to the gluonic flux-tube length of the QQ system. 
Lattice QCD also reproduces this functional form of the 
QQ potential 0, [f|. 

In addition, the lattice QCD calculations reveal that 
this picture holds in three-quark (3Q) systems and multi- 
quark systems. The 3Q potential is obtained from 
quenched lattice QCD as 



V 3 Q(ri,f 2 ,f 3 ) = a 3Q L min ~ ^2 p~r^ 



3Q) 



(2) 



where (r*i, r*2, ^3) arc the coordinates of three quarks and 
L m in is the flux-tube length 0]. The 3Q flux tube has 



the geometry minimally connecting the three quarks, and 
forms the Y-type structure 0, H, 0, E3] • The similar po- 
tentials are obtained in multi-quark systems In the 
flux-tube picture of 3Q and multi-quark systems, the con- 
fining force is a many-body force, reflecting the compli- 
cated gluonic dynamics based on the SU(3) gauge sym- 
metry. The strength of the confinement by these flux 
tubes, i.e., the string tension, is the universal value in 
hadrons, about 0.89 GeV/fm. 

These inter-quark interactions are mainly constructed 
from the gluon dynamics. In particular, as most of pre- 
vious lattice works, the static and quenched calculation 
strictly gives only the gluonic potential. However, since 
not only gluons but also quarks exist in hadrons, the re- 
alistic inter-quark potential would include quark effects. 
The quark effects can be categorized into two types. 

One is the unquenched effect, or the sea quark effect. 
As already known, the sea quark causes "string break- 
ing" and flattens the potential slope in long range. This 
phenomenon is convinced also in lattice QCD [12j | . 

The other is the finite-mass valence quark effect. An 
example of this effect is the relativistic correction to 
the Coulomb interaction, i.e., the Fermi-Breit interac- 
tion [l3j]. As well as the Coulomb potential, we can also 
expect the finite-mass quark effects on the confinement 
potential. In 3Q or multi-quark systems, there would ex- 
ist nontrivial effects reflecting the characteristic flux-tube 
structure. In this paper, we investigate this "motional" 
effect of finite-mass valence quarks. 

To investigate such quark motional effects, we define 
the heavy-heavy- light quark (QQq) potential [l4|. Con- 
sider the QQq system which is constructed with two 
static quarks and one finite-mass quark. The QQq po- 
tential VQQ q (R) is defined as the energy of the QQq sys- 
tem in terms of the inter-heavy-quark distance R. It 
is the inter-two-quark potential in baryons which effec- 
tively includes the light-quark effects. We calculate the 
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FIG. 1: (a)The QQg Wilson loop. The wavy line represents the light-quark propagator and the straight line the heavy-quark 
trajectory. (b)The "wall-to-wall QQq Wilson loop." The gray wavy lines represent the wall-to-wall quark propagator, which is 
the average of all propagators at a fixed time separation. 



QQq potential in lattice QCD, and study the nontrivial 
light-quark effect by comparing it with the static QQ or 
3Q potential. 

As an example of experimental QQq baryons, the dou- 
bly charmed baryon is recently observed by the SELEX 
Collaboration at Fermilab. In 2002, they observed the 
doubly charmed baryon S+ (dec) through a decay process 
S+ — > A+K~tt + They also confirmed another de- 

cay process H+ — ► pD + K~ [16| . However, there is also 
the negative result, for example, by the BABAR Col- 
laboration [r7| . The experimental discovery of the dou- 
bly charmed baryons is still under debate, and more ex- 
perimental information is desired. The doubly charmed 
baryons are also theoretically investigated in lattice QCD 
[IH and other approaches [H, QjJ [2(| ■ 

In this paper, we calculate the QQq potential in 
SU(3) lattice QCD, and investigate the finite-mass va- 
lence quark effect in the inter-quark potential in baryons. 
In Sec. II, we formulate the potential calculation using 
the QQq Wilson loop. In Sec. Ill, we introduce the 
Coulomb gauge fixing for the statistical error reduction, 
and consider its effect on the inter-quark potential. In 
Sec. IV, we show the lattice QCD data and the resulting 
QQq potential, and discuss the finite-mass valence quark 
effect on the QQq potential. Section V is devoted to the 
summary. 



II. LATTICE QCD FORMALISM 

A. QQq Wilson loop 

As the basic potential calculations in lattice QCD, we 
define the QQq Wilson loop for the QQq potential cal- 
culation. As shown in Fig.QJa), the QQq Wilson loop is 
constructed from two "staples" of static-quark trajecto- 
ries and one light-quark propagator. The two staples U 1 
and U u are the path-ordered product of link variables, as 



the static inter-quark potential calculation. The major 
difference from the usual potential calculation is that the 
QQq Wilson loop includes the quark propagator if -1 , 
which represents the light quark moving around. They 
are written as 



1 



W QQq (R,T) = -,e abc e def U l ad Ul l e K-^ 



3! 

U k = Pe l9 k dx " A » 



(k = I, II), 



Kb 



DqDqq a q b e^ K \ 



(3) 

(4) 
(5) 



where the subscripts a, 6, / are color indices. When T 
is large enough, the expectation value of the QQq Wil- 
son loop should depend only on the spatial size R and 
the temporal size T, not on the position of the source 
and sink junction points (i?i and Ru in Fig. [TJa)). The 
resulting QQq potential is written with one parameter i?, 
which is the distance between the two heavy quarks. 

The QQq potential is obtained from the expectation 
value of the QQq Wilson loop as 



,(R) = - lim hn(W QQq (R,T)). 

1 — »-oo 1 



(6) 



The symbol ( ) means the expectation value integrated 
over the gauge field. In practical analysis, the expec- 
tation value is fitted with a single exponential form 
Ce~ VQQqT in large but finite range of T. To estimate 
the suitable fit range of T, the effective mass is defined 
as 



v(R,T) = In 



(W QQq (R,T)) 
(W QQg (R,T + l)Y 



(7) 



If the state is dominated by a single component of the 
potentials, the effective mass is independent of T. The 
plateau region of the effective mass is the indication to 
determine the fit range of T. We calculate the expec- 
tation value in lattice QCD for several values of R, and 
explore a suitable functional form of Vqq 9 (R) . 
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TABLE I: Simulation parameters. The list shows (5 = 2N c /g 2 
and the corresponding lattice spacing a, the sweep numbers 



(JVti 



, N S e P ) of the thermalization and separation for updat- 



ing the gauge fields, the smearing parameters (a, 7V smr ), and 
the clover coefficient c. 



p 


a [fm] 


lattice size 




iVsep 


a 


fVsmr C 


6.0 


0.10 


16 4 


10000 


500 


2.3 


40 1.479 



TABLE II: The correspondence between k and the meson 
masses. The list shows the used gauge configuration num- 
ber N con f, the pion mass m^, the p meson mass m p , and the 
approximate constituent quark mass M q ~ m p /2. The statis- 
tical error is estimated with the jackknife method. 



K 


Nam! 




m p 


Mq 


0.1200 


1000 


1.446(1) 


1.472(2) 


1.5 GeV 


0.1300 


300 


0.900(2) 


0.949(1) 


1 GeV 


0.1340 


300 


0.643(1) 


0.716(1) 


700 MeV 


0.1380 


1000 


0.304(1) 


0.467(2) 


500 MeV 



In Eq. d3J) and Fig. Ufa), the QQq Wilson loop is de- 
fined as a single gauge-invariant loop. For the reduc- 
tion of the statistical error, we actually use the "wall-to- 
wall QQq Wilson loop" . We define the wall-to-wall quark 
propagator as 



K. 



1 



wall V 



(T) = yz K 1 ( n src,n aiDk ,T), 



(8) 



where n src and n s i n k are the spatial sites of the source and 
sink respectively, and V is the number of the spatial lat- 
tice sites. This wall-to-wall propagator is the averaged 
propagator from the whole space at one time to that 
at another time. (The "wall" means the average over 
all spatial sites.) The wall-to-wall QQq Wilson loop is 
constructed by replacing the single quark propagator in 
the QQq Wilson loop with this wall-to-wall quark prop- 
agator. The schematic figure is depicted in Fig. H{b). 
Because such a propagator is independent of the spatial 
position, we can easily sum up the parallel translated 
wall-to-wall QQq Wilson loops in the whole space. This 
summing up drastically suppresses the statistical error 
owing to the large statistics. Without gauge fixing, the 
disconnected components in the wall-to-wall QQq Wil- 
son loops automatically vanish by Elitzur's theorem and 
only the gauge-invariant components remain, and there- 
fore the resulting QQq potential is gauge invariant. 



B. Simulation details 

We perform SU(3) lattice QCD calculation at the 
quenched level. Simulation parameters are summarized 
in Table fl] The gauge action is the standard isotropic 
plaquette action and j3 = 6.0. The corresponding lattice 
spacing a ~ 0.10 fm is determined so as to reproduce 



the string tension <tqq of the QQ potential to be 0.89 
GeV/fm. We use this lattice unit for most part of the 
paper. Lattice volume is 16 4 , and periodic boundary con- 
ditions are imposed on the space-time boundaries. For 
generating the gauge configurations, we use the pseudo- 
heat-bath algorithm and take 10000 sweeps for the ther- 
malization, and 500 sweeps for the separation of each 
configuration. 

For the ground-state component enhancement, we ap- 
ply the smearing method [3, [21| to the spatial link vari- 
ables of the QQq Wilson loop. In the smearing method 
for SU(3) link variables, we iteratively replace the link 
variable U l {n) (i = 1,2,3) with U 2 (n)(e SU(3)) which 
maximizes 



ReTr 



U} (n) { aUi (n) + £ j¥i (U 3 (n) U t (n + j)u] (n + i) 
+U}{n - j)Ui(n - j)U 3 (n + i - . (9) 



Such a replacement does not change the gauge trans- 
formation property of the link variable, and this smear- 
ing method is a gauge-invariant manner. Physically, the 
gauge-invariant smearing method changes a stringy link 
to a spatially-extended flux tube. It does not change 
the physical content such as the potential, and enhances 
the ground-state component so that the statistical error 
is suppressed. The method has two parameters, a real 
parameter a and the iteration number N smr of the re- 
placements, and our choice of a and iV smr is based on 
the static 3Q case @- 

As for the fermion action for the light-quark propaga- 
tor, we adopt the 0(a)-improved Wilson fermion action, 
i.e., the clover fermion action [22j 

•Squark ^ Qn^-nmQm^ (10) 
K nm = $n,m ~ K /,{(1 ~ 7^)^/t(")^n+A,m 

+(1 + 1i*)Ul{n - £)<5n-A,m} 

(11) 

where n and m are the space-time site indices, and other 
indices are omitted. The clover coefficient c in this action 
is given from the mean field value uq of the link variable 
for the tadpole improvement. We determine c and uq 
from the plaquette value P M „ (n) as 



Uq \3 



(12) 



Here ( } means the average of all plaquettes in the gauge 
configurations as the ensemble average. The measured 
mean field value u is 0.87779(2) in our gauge configu- 
rations. To investigate the light-quark-mass dependence, 
we take different four light-quark hopping parameters, 
k = 0.1200, 0.1300, 0.1340, and 0.1380. In Table HI we 



list the results of the meson correlator calculations using 
these hopping parameters. The constituent quark mass 
M q is roughly estimated with the half of the p meson 
mass m p . Our calculations cover the mass region of 0.5 
GeV < M q < 1.5 GeV. 

The Monte Carlo simulations are performed on NEC 
SX-8R at Osaka University. 

III. COULOMB GAUGE FIXING 
A. Statistical error reduction 

In principle, the QQq potential can be calculated with 
the formalism which is mentioned in the Section II. How- 
ever, since the statistical error is severely large, some 
improvements are necessary for statistical error reduc- 
tion. We have tried several methods for the error reduc- 
tion, and adopted the Coulomb gauge fixing procedure. 
In this procedure, the light-quark propagator and the 
heavy-quark trajectories in the wall-to-wall QQq Wilson 
loop is calculated after the Coulomb gauge fixing. The 
Coulomb gauge fixing procedure surprisingly suppresses 
the statistical error in this study. 

Since the wall-to-wall propagator is a spatial average 
of all propagators, the wall-to-wall QQq Wilson loop in- 
cludes disconnected loops, where the two staples and the 
light-quark propagator are not connected at the source 
and sink junction points. The expectation values of such 
disconnected loops vanish in gauge-invariant formalism, 
but do not vanish with the Coulomb gauge. Such con- 
tributions are gauge- variant artifacts by the nonlocal na- 
ture of the Coulomb gauge, and can give the different 
result from the gauge-invariant potential. However, em- 
pirically, the gauge-variant components at the Coulomb 
gauge rapidly decrease and do not affect the resulting 
potential [231 ]. 

In this section, we show several numerical results of QQ 
and 3Q potentials with the Coulomb gauge, and check 
that its influence is small enough for inter-quark poten- 
tial calculations. Other statistical improvement methods 
which we have tried are added to Appendix A. 

B. QQ potential with the Coulomb gauge 

First, we consider the QQ potential with the Coulomb 
gauge [23[ . We set the correlator of two Wilson lines 

W^\R,T) ee ±Tr{L(n,T)L\? 2 ,T)}, (13) 

where L(f, T) is a Wilson line with temporal length T 
and R = \r\ — The expectation value of this correla- 
tor is nonzero with the Coulomb gauge. The QQ poten- 
tial Vqq u1 with the Coulomb gauge is extracted by fitting 
with (W^ ul ) = Cexp(-VCoui T ). The best _ fit regult of 

the on-axis data with V^ ul (R) = crg% ul i? - A%°g l /R + 
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(A) (B) 




T 



FIG. 2: The correctors W^f (left) and W^jf (right) for the 
3Q potential with the Coulomb gauge. R\ and Ru are the 
same in Fig. H^a), an d Rm is the distance between the two 
staples and the Wilson line. 



Cq'q 1 is listed in Table Hill Although the potential seems 
to include some gauge-dependent contributions in small 
T, when the fit range of T is large enough, gauge ar- 
tifacts dump and the potential approaches the physical 
QQ potential. Compared to the gauge-invariant result in 
SU(3) lattice QCD with f3 = 6.0 [6j, the string tension 
a nK l 1S cl° se t° or slightly higher than the physical value 
a Q Q = 0.0534(18), and the Coulomb coefficient A^ 1 is 
almost the same as the physical value Aqq = 0.267(6). 



TABLE III: The best fit parameters of the on-axis QQ poten- 
tial with the Coulomb gauge. The list has different fit ranges 
of T for fitting with {Wq"q 1 ) = Cexp(-Vj| ul T). A dof is the 
degree of freedom. 



fit range 


Coul 


A Coul 

A QQ 


^Coul 

°Q<5 


xV^dof 


T=[l,8] 
T=[2,8] 
T=[3,8] 
T=[4,8] 


0.062(3) 
0.060(2) 
0.058(2) 
0.059(1) 


0.258(10) 
0.254(7) 
0.256(7) 
0.253(5) 


0.609(1) 
0.606(9) 
0.609(2) 
0.605(3) 


3.13 
2.86 
2.69 
0.41 


TABLE IV: The 3Q potential with the Coulomb gauge ex- 
tracted from W^qb'. The best-fit parameters are of the func- 
tion (|16[) with different fit ranges of T. 


fit range 


Coul 
a 3Q 


A Coul 

A 3Q 


/--fCoul 


X 2 /N dot 


T = [l,8] 
T= [2,8] 
T= [3,8] 
T = [4, 8 


0.0498(3) 
0.0474(3) 
0.0466(7) 
0.0482(15) 


0.138(1) 
0.138(1) 
0.141(2) 
0.136(3) 


0.957(3) 
0.959(3) 
0.968(5) 
0.950(11) 


6.14 
8.34 
3.91 
3.95 
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FIG. 3: The effective mass of the waii-to-waii QQq Wilson loop with the Coulomb gauge and the efficiency of the smearing 
method. The lightest quark case, k = 0.1380, is shown. The left graph is the small loop case with R = 2 and (Ri, Ru) = (1, 1) 
and the right is the large loop case with R = 8 and (Ri, Ru) = (4, 4). In each graph, two kinds of data are the result without 
the smearing method (circle) and with the smearing method of (a,JV sml ) = (2.3,40) (square). Thus, the smearing method is 
found to be effective, especially for large R. 



C. 3Q potential with the Coulomb gauge 

Next, we turn to the 3Q case, which is more relevant 
for our QQq potential calculation. As for the 3Q poten- 
tial with the Coulomb gauge, we consider two types of 
correlators. One is the simple correlator of three Wilson 
lines 

WZqa = ie a6c e de/ L ad (fi,T) J L 6e (f 2! T)L c/ (f3,T),(14) 

and the other is the correlator of one Wilson line and two 
"staples" U l and U u 

Wf Q °B = ^e abc e def U l ad UllL cf (r 3 ,T), (15) 

which is a closer geometry to the QQq Wilson loop. The 
schematic figures are shown in Fig. [5] There is not an 
essential difference between the results from these corre- 
lators, but W^Qg 1 has an advantage in accuracy since the 
smearing method is available. 

The numerical results extracted from W^qb' is shown 
in Table [TV] The total number of geometries of W^qq is 
96: < Ri < 4, < Ru < 4, and 1 < R m < 4, except 
for Ri = Ru = 0. The 3Q potential with the Coulomb 
gauge can be written with the same form as the physical 
3Q potential, 

^Coul 

T/Coul _Coul j \ * 3Q ^Coul t-\a\ 

V 3Q - a 3Q ^min ~ 2^ 7^ + °3Q ■ l 10 J 

The 3Q string tension <t^q u1 and the Coulomb coeffi- 
cient A^q u1 are almost the same as the physical values 
a 3Q = 0.0460(4) and A 3Q = 0.1366(11) 0. Considering 
artifacts from the geometrical asymmetry, we have also 
tried a more general fit function 

/iCoul 

T/Coul _ ^Coul T V L /^Coul f-i 7 \ 

V 3Q - a 3Q ^min ~ 2^ |- _ + °3Q > I 1 ') 



where there are five fit parameters. In this case, we have 
found that A^ ul ~ Af 3 oul ~ Ag 5 " 1 and the result is un- 
changed from the fit function (fTB)) . 

From above results, the QQ and 3Q potentials with 
the Coulomb gauge are found to approach the physical 
potentials if the fit range of T is large enough. This would 
hold both in the long-range physics, such as the string 
tension, and in the short-range physics, such as the one- 
gluon-exchange Coulomb potential. Therefore, we can 
expect that the physical QQq potential is approximately 
obtained from the wall-to-wall QQq Wilson loop with the 
Coulomb gauge in the whole region. 



IV. LATTICE QCD RESULTS 
A. QQq potential 

In Fig. [31 we plot typical examples of the effective mass 
of the wall-to-wall QQq Wilson loop with the Coulomb 
gauge. We show the smallest-quark-mass case, n — 
0.1380, where the statistical fluctuation is the largest. In 
the figure, we compare two kinds of data with and with- 
out the gauge-invariant smearing method. The smearing 
method enhances the ground state component of QQq 
potential, especially in the large loop case. Then the ef- 
fective mass is almost flat in T > 3 and thus the ground- 
state component dominates. 

In the region where the effective mass is flat, 
we fit (WQQ q (R,T)) with a single exponential form 
Cexp(— VqqqT), and obtain the QQq potential VQQ q (R). 
The resulting values with k = 0.1380 are listed in Table 
El and k = 0.1200 in Table ED As mentioned above, 
VQQq is almost independent of the fit range of T, and 
then the ground state component dominates. It is also 
confirmed that Vqq 9 does not depend on the position of 
the junction points Ri and Ru separately, but depends 
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k=0.1200 

k=0.1300 o 

k=0.1340 a 

k=0.1380 v 



4 5 
R 

FIG. 4: The lattice QCD data of QQq potential Vqq, with 
the Coulomb gauge. The results of different four hopping 
parameters n are shown. The solid curves are the best-fit 
functions of Eq. (|18p . All the scales are measured in lattice 
unit. 



only on R. Then we can describe the QQq potential as a 
function of the inter-heavy-quark distance R. The same 
arguments hold also in other two n cases. 

For a functional form of VQQ q (R), we consider the 
Coulomb plus linear potential, 



,(R) = a e gR 



Acs 
R 



TABLE V: The lattice QCD results for the QQq potential 
Vqq q with the Coulomb gauge at k = 0.1380. R and (Ri,Rn) 
denote the loop size denned in Fig. [T] The results with differ- 
ent fit ranges of T are also shown. All the values are in lattice 
unit, and the statistical error is estimated with the jackknife 
method. 



a 



cff, 



(18) 



R 


(Ri, Rn) 


T = [4, 8] 


VQQq 

T= [5,8] 


1 


(0,1) 


0.877(2) 


0.873(2) 


2 


(0,2) 


0.971(7) 


0.959(9) 




(1,1) 


0.969(8) 


0.958(10) 


3 


(0,3) 


1.047(4) 


1.045(7) 




(1,2) 


1.045(4) 


1.043(8) 


4 


(0,4) 


1.083(11) 


1.067(17) 




(1,3) 


1.079(10) 


1.063(16) 




(2,2) 


1.078(10) 


1.063(15) 


5 


(0,5) 


1.136(6) 


1.122(3) 




(1,4) 


1.131(6) 


1.117(4) 




(2,3) 


1.130(6) 


1.116(5) 


6 


(0,6) 


1.170(13) 


1.151(24) 




(2,4) 


1.157(16) 


1.136(30) 




(3,3) 


1.157(16) 


1.136(31) 


7 


(0,7) 


1.219(21) 


1.220(50) 




(3,4) 


1.207(24) 


1.209(60) 


8 


(0,8) 


1.262(11) 


1.283(21) 




(4,4) 


1.255(6) 


1.271(10) 



as the analogy of the QQ potential (fTJ) . The subscript 
"efF' means these values effectively including the light- 
quark effect. This simple function is surprisingly suitable 
for Vqq q (R), and the best-fit parameters and the result- 
ing potential form are shown in Table IVIII and Fig. 2J 
respectively. For comparison, the string tension and the 
Coulomb coefficient in the static 3Q potential @ are 

<j 3Q ~ 0.045, ^ 3 q^0.13 (19) 

in the lattice unit at (3 — 6.0 Q- The Coulomb coefficient 
A e g is almost the same value as A^q , which is consistent 
with the short-distance behavior of the 3Q potential in 
perturbative QCD. In contrast, a e g is about 10-20% re- 
duced compared to a^q at n — 0.1300, 1340, and 1380, 
as 

o- c ff < o"3Q- (20) 

In the heaviest 0.1200, the effective string ten- 

sion approximately equals to the string tension. Let us 
call this parameter a e s the "effective string tension" . We 
have found that the effective string tension is smaller 
than the string tension of the static 3Q system. The ef- 
fective string tension strongly depends on k, or the light 
quark mass. 



B. Effective string tension 

In the ground state inter-quark potentials in hadrons, 
the confinement potential is the linear function of the 
flux-tube length L min , and the flux tube forms the shape 
minimally connecting the quarks. As depicted in Fig. 
in the ground state of QQ systems, the flux-tube length 
Lmi n equals to the distance R between the quark and the 
antiquark. The confinement part of the QQ potential 
can be written as a linear function of the inter-quark dis- 
tance i?, and the string tension <Jqq is its proportionality 
coefficient. 

In contrast, in 3Q systems or multi-quark systems, the 
flux tube length L m i n and the inter-quark distances do 
not coincide, and its relation is determined by nontriv- 
ial dynamics of QCD. Thus, the confinement potential is 
a linear function of L m i n but a complicated function of 
the inter-quark distances. In addition, the QQq poten- 
tial effectively includes the heavy-light Coulomb poten- 
tial and the light-quark kinetic energy. Therefore, the 
independence of the QQq potential itself, for example 
that the QQq confinement potential is linear with R, is 
a nontrivial result. 

The lattice QCD results suggest that the QQq confine- 
ment potential is written as the familiar linear potential 
form, but the effective string tension cr G ff is smaller than 



TABLE VI: The lattice QCD results for the QQq potential 
Vqq q with the Coulomb gauge at ft — 0.1200. The notations 
are the same as Table IVl 



R 


(Rl,Rn) 


V QQq 


v QQq 






T= [4,8] 


T=[5,8] 


1 


(0,1) 


1.410(7) 


1.398(4) 


2 


(0,2) 


1.512(10) 


1.492(6) 




(1,1) 


1.510(10) 


1.491(7) 


3 


(0,3) 


1.593(8) 


1.579(8) 




(1,2) 


1.590(7) 


1.577(8) 


4 


(0,4) 


1.637(10) 


1.619(10) 




(1,3) 


1.633(9) 


1.624(8) 




(2,2) 


1.632(9) 


1.614(7) 


5 


(0,5) 


1.694(10) 


1.671(6) 




(1,4) 


1.689(10) 


1.667(5) 




(2,3) 


1.688(9) 


1.667(3) 


6 


(0,6) 


1.724(15) 


1.689(11) 




(2,4) 


1.712(16) 


1.678(16) 




(3,3) 


1.712(16) 


1.678(17) 


7 


(0,7) 


1.782(11) 


1.756(10) 




(3,4) 


1.766(14) 


1.741(23) 


8 


(0,8) 


1.843(8) 


1.846(19) 




(4,4) 


1.834(4) 


1.835(10) 



the static 3Q string tension 03Q . The string tension is the 
proportionality coefficient of the flux-tube length in con- 
finement potentials, and characterizes the confining force 
by the flux tube. The effective string tension is the pro- 
portionality coefficient of the inter-two-quark distance, 
and characterizes the confining force between two quarks 
in hadrons. The deviation between the string tension and 
the effective string tension is considered to originate from 
such a difference, i.e., the geometrical difference between 
the flux-tube length and the inter-quark distance. 

As mentioned before, the functional form of the QQq 
potential is generally nontrivial. Let us consider the more 
detail about its functional form, fn the large R limit, the 
light-quark spreading vertical to the inter-heavy-quark 
direction is negligible compared with R, and the flux- 
tube length is approximately equal to R. This intuitive 



TABLE VII: The best-fit values of a cf f, A g, and C c s in 
Eq. (|18[) . The list also shows the used gauge configuration 
number 7V con f and their % 2 over the degree of freedom TVdof- 



ft 


Aconf 


0" c ff 


A eB 


Ccft 


X 2 /N dol 


0.1200 


1000 


0.045(2) 


0.12(2) 


1.49(2) 


1.31 


0.1300 


300 


0.038(4) 


0.13(2) 


1.23(3) 


1.18 


0.1340 


300 


0.037(4) 


0.13(2) 


1.12(2) 


1.11 


0.1380 


1000 


0.037(2) 


0.13(1) 


0.97(1) 


1.16 




R 



FIG. 5: The schematic figure of the flux-tube length L m m 
and the inter-two-quark distance R. These are equal, i.e., 
L m in = R, in the QQ system (left), and these are not equal, 
i.e., L m i n 7^ R, in the 3Q or QQq system (right). 



consideration suggests that cr e ff would approach to 03Q 
at much larger R, namely the effective string tension is 
some function of R, not a constant. This behavior is 
rather natural, since the QQq confinement potential is 
not necessarily a linear function of R under the nontrivial 
relation between the flux tube length and the inter-quark 
distance. If the QQq confinement potential Vqq^(R) is 
a general function of R, the effective string tension is 
defined as its derivative, 



<Jeff(R) = 



ov QQq 



aR) 



dR 



(21) 



and can depend on R. In the potential model study of 
Ref. (24|, by calculating up to R = 2.4 fm, it is confirmed 
that the effective string tension slightly depends on R. 
The same behavior will be confirmed also in lattice QCD 
with a larger-volume calculation. 

Next we consider the light-quark-mass dependence of 
the QQq potential. When the light-quark mass is larger, 
the spatial extension is more compact and the flux-tube 
length is closer to the inter-heavy-quark distance. In the 
infinite mass limit, the QQq system corresponds to the 
static 3Q system, and the effective string tension equals 
to the string tension. Then the effective string tension 
is an increasing function of the light-quark mass, and 
approaches asymptotically to ct^q in the infinite mass 
limit. We can confirm these behaviors in Table IVlTl 

In Ref. [24[ , the same QQq potential is investigated in a 
non-relativistic potential model, or a quark model. This 
potential model reproduces the present lattice QCD re- 
sult under the same condition of the quark mass and the 
range of R. The potential model can calculate the light- 
quark wave function and the expectation value of the 
QQq flux-tube length. It enables us to understand the 
reduction mechanism of the effective string tension. By 
investigating the relation between the flux-tube length 
and the inter-quark distance R quantitatively, we con- 
firm that a geometrical difference between these is essen- 
tial for the reduction of the effective string tension, as 
conjectured above. 

We have found that, in QQq systems, the effect of the 
finite-mass valence quark reduces the effective string ten- 
sion between the two heavy quarks from the string ten- 



sion of static 3Q systems. This reduction originates with 
the fact that the inter-quark distance differs from the 
flux-tube length in 3Q systems. This is a simple and gen- 
eral property. Our calculation is performed with QQq 
systems for simplicity, however, this simple argument 
would also hold for ordinary baryons, which include three 
finite-mass quarks. Although the finite-mass correction 
is more complicated, the effective string tension can be 
reduced in ordinary baryons, such as a nucleon. Further- 
more, this can be also applied to the multi-quark system 
including light quarks [251 ]. In multi-quark systems, the 
inter-two-quark potential receives the more complicated 
effects of other valence quarks, and the effective string 
tension between the two quarks would be changed from 
the string tension. 





FIG. 6: Examples of the junction choices of the QQq Wilson 
loop. The wavy line represents the light-quark propagator 
and the straight line the heavy-quark trajectory. 



V. SUMMARY 

In summary, we have studied the QQq potential in 
SU(3) lattice QCD, and investigated the role of the finite- 
mass valence quark in the inter-quark potential. For the 
error reduction, we have adopted the Coulomb gauge fix- 
ing and the wall-to-wall quark propagator. From the QQ 
and 3Q potentials with the Coulomb gauge, the Coulomb 
gauge calculation approximately gives the physical po- 
tential in the whole region. 

We have found that the QQq potential is well described 
with a Coulomb plus linear potential, at least in the re- 
gion of R < 0.8 fm. The Coulomb coefficient yl c ff is al- 
most the same as the 3Q case A 3 q, but interestingly, the 
effective string tension a e g is 10-20% reduced from the 
string tension a 3 Q in the static 3Q case. The light-quark 
mass dependence of the potential is also investigated in 
the range of 0.5 GeV < M q < 1.5 GeV. 

The effective string tension is the confining force be- 
tween two heavy quarks in QQq systems. The reduc- 
tion of the effective string tension means that the inter- 
two-quark confining force appears to be weakened by the 
motional effect of the other finite-mass valence quark. 
It originates from the difference between the flux-tube 
length and the inter-quark distance, and reflects the char- 
acteristic flux-tube structure of baryons. 

This reduction of the inter-two-quark confinement 
force is conjectured to be a general property not only 
for QQq systems but also for ordinary baryons. Also 
in multi-quark hadrons, we can expect similar or more 
complicated effects on the inter-two-quark potential by 
the finite-mass valence quark. The quark confinement 
is a fundamental property for hadrons, and its change 
would be important for broad fields relating quark- 
hadron physics. 
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APPENDIX A: OTHER TRIALS ON 
STATISTICAL IMPROVEMENT 

We have tried several statistical improvement tech- 
niques for calculating the QQq potential. These tech- 
niques are useful for error reduction to some extent, but 
not enough to calculate the QQq potential in gauge- 
invariant manner. The techniques are not used for the 
final result of the Coulomb gauge calculation. We intro- 
duce them briefly in this Appendix. 



1. Multi-hit procedure 

We have applied the multi-hit procedure [III [27| • The 
temporal link variables are replaced with the mean-field 
value of the neighboring link variables. It is realized by 
replacing the temporal link variable U4 (n) in the staples 
with C/4(n)(G SU(3)) which maximizes 

Ej ReTr \^l{n){Uj{n)Ui{n + j)Uj(n + 4) 

+U](n - j)U 4 (n - j)Uj{n + 4 - j)}] • (Al) 

2. Average of the junction points 
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In Fig. [TJ the positions (Ri,Ru) of the junction points 
are the same at the source and sink of the QQq Wilson 
loop. We can take (Ri, Ru) at the source and sink inde- 
pendently, and average all the combinations of (Ri, Ru) 
with fixed R. This improvement increases the statistics 
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R 



FIG. 7: The QQq potential Vqq 9 of the different gauge 
choices: (A) gauge invariant with k = 0.1200 and JV CO nf = 
4000, (B) Coulomb gauge with k = 0.1200 and AW = 1000, 
(C) Coulomb gauge with « = 0.1380 and N cont = 1000, (D) 
Landau gauge with k = 0.1380 and Ar con f = 200. The solid 
curves are the best-fit functions of Eq. (|18|) . 



TABLE VIII: The results of the different gauge choice. The 
list shows the configuration number N con f, and the best-fit 
values of a e s and A e s in Eq. (|18[) . A statistical improvement 
is adopted at the gauge-invariant calculation. 



gauge 


K 


N r 
1 * conr 


Cett 


A eB 


gauge inv. 


0.1200 


4000 


0.051(3) 


0.15(1) 


Coulomb 


0.1200 


1000 


0.045(2) 


0.12(2) 


Coulomb 


0.1380 


1000 


0.037(2) 


0.13(1) 


Landau 


0.1380 


200 


0.039(2) 


0.13(2) 



by (R + l) 2 times, and effective in large R, where the 
statistical error is severe. 
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3. More average of the junction points 

The positions of the junction points can be taken more 
arbitrarily. For example, as depicted in Fig. [5J we have 
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the staples, and averaged such contributions with fixed 
R. In addition, we would be able to take more arbitrary 
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APPENDIX B: GAUGE DEPENDENCE 

We show here the results of the QQq potential without 
gauge fixing or with another gauge. 
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ror is relatively small and the computing time for the 
quark propagator is fairly short. For further statistical 
improvement, the junction-average procedure explained 
in Appendix A2 is adopted for the gauge-invariant cal- 
culation. Nonetheless, the statistical error is still large, 
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fm. The result is shown in Fig. [7] and Table IVIIll Note 
that the listed error in Table IVIIll is only the statistical 
error, and that the gauge-invariant calculation involves 
the systematic error from the fit-range dependence. The 
gauge-invariant result seems to be close to the Coulomb 
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the large statistical and systematic error. 
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